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Paramagnonlike excitations and spin diffusion in magnetic resonance studies 

of copper oxide superconductors 

Igor A. Larionov^ 

Magnetic Radiospectroscopy Laboratory, Department of Physics, Kazan State University, 4^0008 Kazan, Russia 

The relaxation function theory for a doped two-dimensional Heisenberg antiferromagnetic system 
in the paramagnetic state for all wave vectors through the Brillouin zone is presented in view of 
low frequency response of high- Tc copper oxide superconductors. We deduced the regions of long 
lifetime [T <400(1 — 4a;) K] and "overdamped" [T >700(1 — 4a;) K] paramagnonlike excitations in the 
temperature (T)-doping index (x) phase diagram from plane oxygen nuclear spin-lattice relaxation 
rate "(1/Ti) data in up to optimally doped La2-i:Sra;Cu04 thus providing the regimes for the spin 
wave concept and the "overdamped" mode. 



I. INTRODUCTION 

Plane copper oxide high-temperature superconductors 
(high- Tc) are the doped 8—1/2 two dimensional Heisen- 
berg antiferromagnetic (2DHAF) systems.^ In the carrier 
free regime, the elementary excitations are spin waves,— 
magnons in the quasiparticle language, a concept widely 
known and thoroughly investigated in the pasti^ There- 
fore it is tempting to consider the doped 2DHAF sys- 
tems in terms of magnonlike excitations (strictly speak- 
ing, the paramagnon, a notation used for spin fluctuations 
in representation of damped spin waves) and the magnon 
lifetime is characterized by the damping of spin waves. 
Then the questions arise: what happens with spin waves 
when we dope the system? What are the elementary ex- 
citations - are they still paramagnonlike? The motion of 
charge carriers even in the optimally doped (maximum Tc) 
high- Tc is known to take place in the presence of strong 
AF fluctuations^ and spin waves in 2DHAF systems per- 
sist even without long range order in the paramagnetic 
state^i^i^i^ so the questions make sense, but, finally, what 
one can say about the lifetime of these excitations when 
we dope the system? 

From experimental point of view the spin-wavelike fea- 
tures have been revealed^, by neutron scattering (NS) even 
in the nearly optimally doped YBa2Cu3 06.85 and, con- 
trary to the predictions within the weak coupling theory^ 
no isotope effect on the "resonance peak" (RP) frequency 
have been observed^ in YBa2Cu3 06.89- The RP phe- 
nomenon disappears in the overdoped phaseiS, thus raising 
questions about its appearance within the weak coupling 
theoriesii and AF spin excitations disappear in overdoped 
Lai.7Sro.3Cu04 leading to the conclusion that "the AF 
spin correlations in superconducting samples must be ves- 
tiges of the parent insulator" Moreover, NS data— 
in underdoped YBa2Cu3 06.35 with Tc=18 K shows the 
commensurate AF short range order, no well-defined reso- 
nant mode, and similarities between the low-energy mag- 
netic excitations in YBa2Cu3 06.35 and carrier free insu- 
lating 2DHAF, e.g., YBa2Cu306.i5 and La2Cu04J^ The 
spin diffusive contribution ^"^ {l/Ti)oif f to plane oxygen 
nuclear spin-lattice relaxation rate '^^(1/Ti) cannot be ex- 
cluded solely by failure to detecl— the changes in ^''(l/Ti) 
by varying the nuclear magnetic resonance (NMR) fre- 



quency Ld since ^^(l/Ti)oi// varies rather weaklyii with 

LJ. 

In this paper we use the Mori-Zwanzig projection opera- 
tor procedure^^'^^ and thus we are unprejudiced regarding 
the role of q « (spin diffusion) and Q ~ (tt, tt) wave vec- 
tors in the imaginary part of dynamic spin susceptibility 
x"(q, w) of doped 2DIIAF system which is especially im- 
portant for ^^(1/Ti). Spin diffusion is a spatial smoothing 
of heterogeneous spin polarization in a system of localized 
magnetic moments and in the presence of strong damp- 
ing (short magnon lifetime) the spin dynamics changes 
from wavelike to diffusive. The approximations we use 
for the relaxation function are within the Markovian ap- 
proximation and "by itself the Markovian situation can be 
valid even in the absence of any picture of the system in 
terms of well-defined excitations" We will emphasize the 
spin- wavelike features in x" (Q: of copper oxide high- Tc 
and extract the lifetime of spin-wavelike excitations from 
"(1/Ti) data. 



II. BASIC RELATIONS 

We start from the t-J Hamiltonian^S known as the mini- 
mal model for the electronic properties of high- Tc cuprates 

i,j,<y i>j 

written in terms of the Hubbard operators ATf " that cre- 
ate an electron with spin a at site i and are spin- 1/2 
operators. Here, the hopping integral tij=t between the 
nearest neighbors (NN) describes the motion of electrons 
causing a change in their spins and J=0.12 eV is the 
NN AF coupling constant. The spin and density opera- 
tors are defined as follows: S'f ^Xf^, = (1/2)^;^ crX^ , 
ni—'Y^^Xf'^, (ct=— ct), with the standard normalization 
X00+X+++X-- = 1. 

We formulate our study of the spin ffuctuations follow- 
ing Mori,-'^^ who showed it's efficiency for both the clas- 
sical (and essential equivalence to Brownian motion) and 
quantum (e.g., Heisenberg systems of arbitrary dimension) 
many body systems.— The time evolution of a dynamical 
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variable S'^(t), say, is given by the equation of motion, 



dT 



zCSUt)^[H,^j,S^{t)], (2) 



where the Liouville operator C in the quantal case rep- 
resents the commutator with the Hamiltonian. The pro- 
jection of the vector S^{t) onto the = S^{t=0) axis, 
VqS^{t) = 7?.(k, t)5'^, defines the hnear projection Her- 
mitian operator Vq. One may separate S^{t) into the 
projective and vertical components S^{t) = 7?.(k, t)S'^ + 
(1 — 'Pa)S^{T) with respect to the axis, where 
7^(k,r) = (S'^(t), {Sli^)*){S^, ('5'lk)*)"^ is the relaxation 
fimction in the inner-product notation: {S^{t), (5^^)*) = 

ksT J^^'"'"'^ dg{ee"S^{T)e-^"{S\)*), and the angular 
brackets denote the thermal average. 

One may construct a continued fraction representation 
for the Laplace transform of the relaxation function, for 
which Lovesey and Meserve2i2i used a three pole approx- 
imation, n^{k,s)^J^ dr e-^^7^(k,T)«l/{s + Ajjls + 
^2k/('* + I/t'Ic)]}, with a cutoff characteristic time Tk = 
^2/(7rA2k), by arguing that S^{t) fluctuations are weakly 
affected by higher order random forces. For the relaxation 
shape function J^{'k,uj)=Re[R.^{k,iuj)]/n, this gives 



(3) 



where A^^ and arc related to the frequency moments, 
K> = fZodiJiv^J^iKu;) = (l/^")[d"7^(k,r)/dr"],^„, 

of 7^(k,r) as A?k - H), A^ = (H) / H)) ~ 
(wk) for T > Tk. Note that J^(k, w) is normalized to 
unity (iwJF(k, uj) ^ 1 and is even in both k and 
CO. The expression for the second moment, (wk) — 
i{[S^,Sl^,])/xik) = -(8Jci-4t,//ri)(l-7k)/x(k), is 
compact, while (wk) = St]^]) /yik) is cumbersome 

and is not reproduced here (see Ref. [ijfor details). 

The static spin susceptibility has been derived within 
the t-J model in the overall temperature and doping 
ranges^ 



X(k) 



4ki| 



(4) 



and has the same structure as in the isotropic spin- 
wave theory.— The parameter g+ is related to corre- 
lation length ^ via the expression £^/a^l/{2y/g^ — 1), 
where a=3.8 Ais a lattice unit. The transfer 
amplitude between the NN is given by: Ti = 
-(1/4) Ep (^f"^»>p)=PEk7k/^ the index p runs over 
NN, 7k=(l/4) E^exp(ik/9)=(l/2)(cosfca;a-|-cos/cj^a), and 

/k = [exp(— i?k + p)/kBT + 1]^^ is the Fermi function of 
holes. The number of extra holes, due to doping, 6, per 
one plane Cu'^"*' , can be identified with the Sr content x in 
La2_2:Sr2,Cu04. The chemical potential /i is related to S 
by S—pJ2k flij with p={l+5)/2. The excitation spectrum 
of holes is given by E]s_=Atef f^fic, where the hoppings, t, 
are affected by electronic and AF spin-spin correlations 
ci, resulting in effective values,— for which we set 



tgff = 5 J/0. 2 in order to match the insulator-metal tran- 
sition. 

For low temperature behavior we use the expression, re- 
sulting in effective correlation length f , given b y^^i^^'^^ 



eff 



(5) 



Here, ^ is affected by doped holes, in contrast with 
the Keimer et al^ empirical equation, where ^ is given 
by the Hasenfratz-Niedermayer formula and there was 
no infiuence of the hole subsystem on ^. Thus from 
now on we replace ^ by ^e//- For doped systems 
we use the explicit expression^^ for f which is much 
more complicated compared with simple relation ^/a ~ 
{J^/~gZ/kBT)exTp{2TTps/kBT), valid for carrier free or 
lightly doped systemSfi^i^ In the best fit of ^e// to ex- 
perimental dat a^^i^^ the relation i^o = a/n^x is most 
suitediiiii which one may attribute to stripe picture, 
where n^=2 for x <0.05 and nj=l near optimal (x ~0.15) 
doping. The results of the calculations are summarized in 
Table I. We consider here the case of La2-a:Sr3;Cu04 with 
the simplest crystalline structure for brevity and luck to 
thorough experimental data set. 



III. PLANE OXYGEN NUCLEAR SPIN-LATTICE 
RELAXATION 



The plane oxygen nuclear spin-lattice relaxation rate 
^^(1/Ti) has three contributions: 

"(1/Ti) = "il/T,),^,+'\l/Ti)Korr + '\l/Ti)D^ff. (6) 

The contribution from spin-wavelike excitations is given 
by 
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(1/Tl). 



2kBT 

Wo 



''FikfxUKu^o), (7) 



|k|>l/«e// 



where luq = 2ttx52 MHz ~2. 15x10""* meV (< T, J) is the 
measuring NMR frequency at 9 Tesla. The quantization 
axis is along the crystal c axis and the wave vector depen- 
dent hyperfinc form factor for plane *^0 sites is given by 
"F(k)2 = 26*2(1 + 7k), with C = 2.8 X 10"^ eV.-^ 

The contribution from itinerant holes, of Korringa type, 
^"^ {l/Ti)Korr = KkT, should grow with doping x and will 
be the adjustable parameter. The contribution from spin 
diffusion (small wave vectors k < 1/Ce//) rnay be calcu- 
lated from general physical grounds, namely, the linear re- 
sponse theory, hydrodynamics, and fluctuation-dissipation 
theorem^^ (see also Ref. 17) 
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(l/Ti) 



Dtff 



^F{QfkBTa^x{k ^ 0) 



A, 



(8) 



where A=[l/(47r)] ln[l -I- i'Vl^o^l//)] and the calculated 
values of spin diffusion constant, l), are given in Table I. 

Since the relaxation function can be understood within 
the spin-wave framework,— the temperature and doping 
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TABLE I. The calculated in the T ^0 hmit NN AF spin-spin 
correlation function Ci = (l/4)^^{5'j^5j?^p), the parameter g-, 
the spin stiffness constant ps using the expressions and the 
procedure as described in Refs. and [23 . the calculated spin 
diffusion constant, D, following Ref. [l3 together with the 
values of Korringa-type contribution constant Kk and the 
spin-wavelike damping renormalization constant Tr as 
extracted from comparison with ^''(l/Ti) NMR data. 



X 


Cl 


9- 


2-Kps/J D/Ja^ 


Kk, (sK)-i 


r., K- 


3 


0.025 


-0.1133 


4.102 


0.36 


2.60 


0.023 


4.1x10 


-9 


0.035 


-0.1115 


4.060 


0.35 


2.54 


0.024 


5.5x10 


-9 


0.05 


-0.1018 


3.827 


0.285 


2.47 


0.051 


7.5x10 


-9 


0.1f5 


-0.0758 


3.252 


0.2 


3.51 


0.147 


32x10" 


-9 


0.15 


-0.0617 2.947 


0.13 


3.81 


0.215 


41x10" 
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dependence of the damping of the spin-wavelike excita- 
tions may be studied further. The spin-wavelike disper- 
sion, rcnormalized by interactions, is given by the relax- 
ation functionf^ 



2 duj ti>jr(k, cj), 
Jo 



(9) 



where the integration over lo in Eq. ([9]) has been performed 
analytically and exactly. We assume the Lorentzian form 
of the imaginary part of the dynamic spin susceptibility. 



XL (k,c.) 



X (k) t^Fk 



x(k)tjFk 



(10) 



for k around (tt, tt). We accept to the leading order the cu- 
bic temperature dependence^fiiii for the damping of spin- 
wavelike excitations rk=rrT'^77k, where the wave vector 



dependence is given by 77k='\/(ij-'k) ~ ("^k™)^ 



IV. RESULTS AND DISCUSSION 

Figure [1] shows the plane oxygen ^"^{l/Ti) fitted by 
Eq. ([6]) with two adjustable parameters: Kk and F^, 
which values are given in Table I. The quality of the fit 
is very good, which we treat as the validity of our the- 
ory. The importance of ^"^ [l/TijKorr and ^'^{l/Ti)oiff in 
plane oxygen in contras t with plane copper 

is due to the filtering of the (7r,7r) contribution 
by the plane oxygen hyperfine form factor. Obviously, in 
the absence of ^"^ {l/Ti)Korr and ^'^ {1/Ti)£)if f , it is hard 
to explain the measured at large x with any form 

of the damping function. In general, the damping grows 
with doping it should. It should be emphasized that 
the increase of ^^(1/Ti) with temperature is caused by the 
increase of the damping, Fk- The low T region of lightly 
damped paramagnonlike excitations, where the data may 
be explained by the theory that neglects the damping, ■'^^ 
is quantified through the relation TrT^ <0.2. 

Figure [2 shows the temperature (T)-doping index (a:) 
phase diagram with the spin-wavelike damping regimes 
deduced from plane oxygen nuclear spin-lattice relaxation 
rate ^^(l/Ti) data. It is tempting to speculate that the 
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FIG. 1: (Color online) Temperature and doping behavior of 
plane oxygen ^^(1/TiT) for La2-a:Sra,Cu04. Dashed lines show 
the spin diffusive contribution. Thin solid curves are the fits by 
Eq. (|6j to NMR datai^i^ as described in the text. Thick solid 
lines separate schematically the regions of lightly damped (low 
T region) , damped, and overdamped magnon-like excitations. 



doping and temperature behavior of these curves resem- 
bles the characteristic "pseudogap" temperatures. 

The form of x'l (k, ^) in Eq. pH)) gives the commensu- 
rate response at low lu and the incommensurate response 
at high uj in agreement with NS studies in the lightly doped 
regimeii^iii^ Very recently. Stock et ali^'^ reported the 
evidence for spin waves from NS studies of underdoped 
YBa2Cu306.35 with Tc=18 K, where the magnetic excita- 
tions are very similar to that of carrier free 2DHAF sys- 
tems. These observations, together with the undoubtful 
evidence for disappearance of AF spin excitations in the 
overdoped regime," show that AF spin excitations in the 
overall doping range of copper oxide high- Tc emanate from 
those of the parent insulator, i.e., the spin waves. 

The wavevector dependence of spin- wavelike dispersion 
Wk™ and damping 77k is shown in Figs.[3Ua) andlS^b), re- 
spectively, for various doping levels. Both lo^ and ryk 
show negligible temperature dependence below T<J/2 
except the region around (7r,7r). For a;=0, our calcu- 



lated magnon energy, , at (7r,0) is 3% lower than at 
(7r/2,7r/2) in qualitative agreement with Monte Carlo sim- 
ulations and series expansion calculations^ and is similar 
to that in Sr2Cu304Cl2,'^^ however, is a bit different from 
NS data in La2Cu04fi^ The dispersion of remains 
approximately the same and the wave vector dependence 
of the damping function, ?7k, in contrast, possesses sig- 
nificant changes with doping. We emphasize that below 
T w400(l-4x) K, where F^T^ <0.2, the damping Fk is 
much smaller compared with t^'k'^, thus the spin- wavelike 
excitations are indeed well defined (long lifetime). Above 
T «700(l-4a;) K, where F^T^ >1, the damping, Fk, be- 
comes compatible with lo^ , and grows further with T and 
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FIG. 4: (Color online) The calculated field (frequency) depen- 
dence of ^^(l/TiT) (solid curve). The dashed curve shows the 
spin diffusive contribution. The error bar shows the size of the 
symbol (error bar) in Fig. [1] Short solid horizontal lines indi- 
cate the calculated values of ^^(l/TiT) at the fields 9 T and 
14.1 T. 



FIG. 2: (Color online) La2-2,Sra;Cu04 phase diagram. Solid 
curves show the Neel temperature (Tjv), the superconduct- 
ing dome (Tc), and dR/dT=0 curve indicates the gradual 
crossover from insulating to metallic behaviour from resistivity 
measurements.— Lower and upper shaded lines with circles have 
been extracted from ^^(l/TiT) data shown in Fig. [1] with the 
conditions r,.r='=0.2 and F^T^^l, where F^ is given in Table I, 
may be approximated by T f«400(l-4a;) and T f«700(l-4a;), 
respectively, and separate the regions of lightly damped (long 
lifetime, well defined) , damped, and overdamped paramagnon- 
like excitations. The asterisks mark approximately the corre- 
sponding temperatures as extracted from the data analysis^^ 
for Sr2Cu02Cl2. 
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( kx/TT, ky/jl ) 

FIG. 3: (Color online) (a) Spin-wavelike dispersion u)^ and 
(b) damping function rjk at various dopings along (i,i)-(l,l)- 
(l,0)-(i,i)-(0,0)-(l,0) route in the Brillouin zone. 



X thus providing the "overdamped" mode region. '^^'^^ 

It should be mentioned that the so-called anomalous 
"1/ 8" doping problem can be viewed as a consequence of 
the particular parameter set at x ~ 1/8 within an ex- 
tended t-J model with the Coulomb repulsion between 
the NN and taking into account the polarization of NN 
copper spins around copper-oxygen singlet that gives the 
in-phase domain structure for the "stripe" picturci^ The 
in-phase domain is favorable, also in view of the exper- 
imental data,— compared with the anti-phase domain 
model4S This particular case is beyond our present con- 
sideration because of its narrowness. 

Figure [4] and Eq. [8] show that in doped 2DHAF system 
the spin diffusive contribution ( 1 /Ti ) di// scales with the 
NMR frequency to as ln{constJ/Lu), which is very weak in 
view of the extraordinary large superexchange coupling 
constant J. We argue that at low doping level (a;«0.03 
per Cu site in La2_a;Sra;Cu04) and temperature T<300 K, 
i^(l/ri) is strongly affected by "(l/Ti)Di// luck of the 
filtering of (tt, tt) contribution by the oxygen formfactor. 
Only a new accurate NMR experiments at very low and 
very high magnetic fields may uncover ^"^ {1/Ti) oif f con- 
tribution to "(l/Ti). 



V. CONCLUSION 

We applied the relaxation function theory for doped 
2DHAF system in the paramagnetic state and deduced 
the lifetime of spin- wavelike excitations, its evolution with 
doping, from plane oxygen ^^(1/Ti) NMR data in the un- 
derdoped high- Tc layered copper oxides. It is shown that 
the spin-wavelike theory is able to reproduce the main fea- 
tures of low frequency spin dynamics in the normal state 
of high- Tc cuprates as observed experimentally. We iden- 
tified the regions of long lifetime [T <400(1— 4a:) K] and 
"overdamped" [T >700(1— 4.t) K] paramagnonlike excita- 
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tions in the temperature (r)-doping index (x) phase dia- 
gram in up to optimally doped La2-2:Srj,Cu04. The re- 
sults indicate that spin-wavelike excitations are indeed a 
good description of the quasiparticle excitations even for 
strongly doped high- Tc layered cuprates at low tempera- 
tures, T <400(l-4a;) K. 
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